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We analyze the distribution of the chromoelectric field generated by a static quark-antiquark pair 
in the SU(3) vacuum. We find that the transverse profile of the flux tube resembles the dual 
version of the Abrikosov vortex field distribution and give an estimate of the London penetration 
length in the confined vacuum. 
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1. Introduction 



Color confinement in Quantum Chromo-Dynamics (QCD) is a long-distance behavior whose 
understanding continues to be a challenge for theoretical physics [[TJ, 01. Tube-like structures emerge 
by analyzing the chromoelectric field between static quarks [g, § |, |, 0, g, g, 0, 0, 0' 



15, 16, 17, 18, 19]. Such tube-like structures naturally lead to linear potential and consequently to 



a "phenomenological" understanding of color confinement. To explain the formation of chromo- 
electric flux tubes in QCD vacuum, 't Hooft [20] and Mandelstam proposed the hypothesis 
that QCD vacuum behaves like a coherent state of color magnetic monopoles, which leads to the 



picture of QCD vacuum a magnetic (dual) superconductor [ |22| ] . According to this picture the ob- 
served color flux tubes are naturally accounted for by the (dual) Meissner effect, in analogy with 
the formation of Abrikosov tubes in the usual superconductivity [^]. Even if the 't Hooft con- 
struction does not explain the dynamical formation of color magnetic monopoles, many lattice 
calculations [^4|, 25, 26, 27, 28, 2^, pC|, pi, 32] have given numerical evidence in favor of magnetic 
monopole condensation in the QCD vacuum. 

On the other side, magnetic monopole condensation could be the consequence rather than the 



origin of the confinement mechanism [|33[]. Even in this case, however, the dual superconductivity 
picture provides us with a "phenomenological" frame for the analysis of tube-like structure in the 
QCD vacuum. 



The outcome of previous studies [ ]12| , |13J, |14j, |1_5J, |16p of the SU(2) confining vacuum was the 
following: (1) presence in lattice configurations of color flux tubes made up by the chromoelec- 
tric fields directed along the line joining a static quark-antiquark pair; (2) transverse size of the 
chromoelectric flux tube interpreted as the London penetration length in the Meissner effect; (3) 
penetration length measured both in the maximal Abelian gauge and without gauge fixing, with 
compatible results, thus supporting gauge-invariance; (4) determination of the string tension as the 
energy stored into the flux tube per unit length, in good agreement with the results in the literature. 

The aim of the present work is to extend the analysis to the more interesting case of the SU(3) 
gauge theory. 



2. Color fields on the lattice 

We use a connected correlator (Fig. [j])(left) to explore the field configurations produced by a 
static quark-antiquark pair (N is the number of colors) [g, g, f34]> |35|] : 
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The configuration with plaquette parallel to the Wilson loop leads to chromoelectric field longitu- 
dinal to the axis defined by the static quarks. 
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Figure 1: (Left) The connected correlator (2.1) between the plaquette U p and the Wilson loop. The subtrac- 
tion appearing in the definition of correlator is not explicitly drawn. (Right) Longitudinal component of the 
chromoelectric field versus the distance x t at /3 = 5.9 after 10 cooling steps. 



2.1 The case of SU(2) gauge theory 



The case of the SU(2) gauge theory was studied in [ TO , T^, 13, 14, 15, 16]. The main result 
was that the flux tube is almost completely formed by the longitudinal chromoelectric field, Ei, 
which is constant along the flux and rapidly decreasing in the transverse direction x t . In the dual 
Meissner effect inteipretation, the transverse shape of Ei is the dual version of the Abrikosov vortex 
field distribution [Jl(J |l2|, ||, 0. IH> Hi and therefore must obey 



Ei{x t ) = —Ll z K (nx t ) , 
2K 



x t >0 



(2.4) 



where <I> is the external flux and X = 1 j\l is the London penetration length (this is valid for A 3> § , 
with £ the coherence length of a type-II superconductor). 



In the past numerical study (lattices 16 4 , 20 4 and 24 4 , statistics 20-100) [16] approximate 
scaling was found, n/y/o = 4.04(18) (i.e. A = 0.118(5) fm for yfc = 420 MeV). In this work 
(20 4 lattice, statistics 1000), we find n/y/o = 4.21(16) - see Ref. [Q for details. 



2.2 The case of SU(3) gauge theory 

The main motivation for repeating the study in SU(3) is to verify the scaling of /» with the 
string tension and to compare the resulting determination of ll/yfG with SU(2). This result should 
provide us with important reference values, that any approach aiming at explaining confinement 
should be able to accommodate. 

We performed numerical simulations with the Wilson action and periodic boundary conditions, 
using a the Cabibbo-Marinari algorithm [[37|], combined with overrelaxation on SU(2) subgroups. 
The summary of /3 values, lattice size, Wilson loop size and statistics is given in Table |I} The 
lattice size L has been chosen such that the combination L^fo > 4. The size of the Wilson loop 
entering the definition of the operator given in Eq. (2.1) has been fixed at L/2 — 2a. In order to 



reduce the autocorrelation time, measurements were taken after 10 updatings. The error analysis 
was performed by the jackknife method over bins at different blocking levels. 

In order to reduce the quantum fluctuations we adopted the controlled cooling algorithm. It is 



known [38] that by cooling in a smooth way equilibrium configurations, quantum fluctuations are 
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p 


lattice 


Wilson loop 


statistics 


5.90 


18 4 


7x7 


5.k 


6.00 


20 4 


8x8 


4.5k 


6.05 


22 4 


9x9 


3.6k 


6.10 


24 4 


10 x 10 


2.4k 



Table 1: Summary of the Monte Carlo simulations. 



reduced by a few order of magnitude, while the string tension survives and shows a plateau. We 
shall show below that the penetration length behaves in a similar way. The details of the cooling 
procedure are described in Ref. JTB] ] for the case of SU(2). Here we adapted the procedure to the 
case of SU(3), by applying successively this algorithm to various SU(2) subgroups. The control 
parameter 8 was fixed at the value 0.0354, as in Ref. [16]. 



A novelty with respect to the study of Ref. fllq ] is related with the construction of the lattice 
operator given in Eq. (2T). If the Wilson loop lies on the plane, say, 1-2, then the Schwinger line 
can leave the plane 1-2 in the direction, say, 3; before attaching the plaquette to the Schwinger line, 
the latter can be prolongated further in the direction 4, by one or two links. In this way, by varying 
the length of the Schwinger line in the direction 3, one can obtain a large set of distances x t /a 
between the center of the plaquette and the center of the Wilson loop, both integer and non-integer. 
On each configuration we averaged over all possible directions for the relative orientation of the 
Wilson loop to the Schwinger line. 

The general strategy underlying this work is the following: (1) for each j8 we generate an en- 
semble of thermalized configurations and, correspondingly, ensembles of "cooled" configurations 
after a number of cooling steps ranging from 5 to 16; (2) for different values of the distance x t , the 
longitudinal component of the chromoelectric field, averaged over each cooled ensemble of con- 
figurations, is then determined by means of the operator (2A), with the help of Eq. (2.3) (see, for 
example, Fig. fright), which shows Ei{x t ) averaged over the ensemble at /3 = 5.90 after 10 cooling 
steps); (3) for each cooling step, data for Ei(x t ) are fitted with the function given in Eq. Q and 
the parameters /J. and <I> are extracted; (4) a plateau is then searched in the plot for jU and <I> versus 
the cooling step. 

In Table ^ we report the results for a}i of the fit at the four j8 values considered in this work for 
one selected cooling step. The table with the results for the other fit parameter, <I>, can be found in 
Ref. [36]. When the fit is done on all available data for E x (x t ), above a certain x t . m i n , the # 2 /d.o.f. 
is very high, thus reflecting the wiggling of data due to the inclusion of non-integer distances x t /a. 
When the fit is restricted to integer values of x t /a, the # 2 /d.o.f. turns out to be very reasonable. 
Remarkably, the resulting parameters obtained with the two fitting procedures agree very well. 

In Figs. ^ we show the behavior of ajx and <1> with the cooling step at jS = 6.05. Similar figures 
for the other values of j3 are given in Ref. [36]. A short plateau is visible, except for the case of \i 
at p = 5.90. We take as "plateau" value for /I the value corresponding to the number of cooling 
steps given in the second column of Table ||. 

Finally, we studied the scaling of the "plateau" values of ajJ. with the string tension. For this 
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Figure 2: (Left) The inverse of the penetration length aji at j3 = 6.05 versus the cooling step. Data are ob- 



tained by fitting the transverse profile of the longitudinal chromoelectric field with the function (2.4); circles 
correspond to fit to all available data of E[(x t ) starting from a certain x, m i n , while squares correspond to fit 
of E[(x t ) for integer values of x t /a. (Right) The same for the amplitude of the longitudinal chromoelectric 
field 4>. 
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cooling step 


ayi 


£ 2 /d.o.f. 


^-(.min / " 


data set 


5.90 


10 


0.5577(12) 


626. 
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all data 


6.00 


9 


0.51015(92) 


383. 
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all data 


6.05 


10 


0.4730(13) 


133. 
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all data 


6.10 


10 


0.4357(20) 


27. 
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all data 


5.90 


10 


0.5557(40) 


1.22 


7 


integer x t /a 


6.00 
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0.5099(28) 


2.56 


9 


integer x t /a 


6.05 


10 


0.4735(39) 


1.08 


8 


integer x t /a 


6.10 


10 


0.4349(56) 


0.25 


8 


integer x t /a 



Table 2: Summary of the fit values for apt. 



purpose, we have expressed these values of a\i in units of ^/o, using the parameterization 
aV*(g) = fsu(3) (g 2 ) [1 + 0.2731 a 2 (g) - 0.01545a 4 (g) + 0.01975 a 6 (g)]/0.01364 , 



fsu(3)(g 2 ) 



fsi/(3)(g 2 (P 

f sm (g 2 ) = (bog 2 y b>/2bl 



6)) 



exp 



1 



5.6 < j3 < 6.5 



11 



102 



(2.5) 



(2.6) 



2b g 2 J ' u (4tt) 2 ' 1 (4tt) 4 ' 
given in Ref. @. 

Figure || suggests that the ratio /I / yfo displays a nice plateau in jS, as soon as /3 is larger than 
6. The scaling of pL is a natural consequence of the fact that the penetration length is a physical 



quantity related to the size D of the flux tube [10, |2|], D ~ 2/jU. We get as estimate for the 
penetration length in SU(3) gauge theory, yi / *fa = 2.325(5), corresponding to jJ. = 0.977(2) GeV. 
We observe that this value is in nice agreement with the determinations of Ref. [ffij, obtained by 
using correlators of plaquette and Wilson loops not connected by the Schwinger line, thus leading 
to the (more noisy) squared chromoelectric and chromomagnetic fields. 
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Figure 3: Scaling of the inverse London penetration length with \fa versus j3. Data have been slightly 
shifted on the horizontal axis for the sake of readability. 



We note that the ratio between the penetration lengths for the SU(2) gauge theory and the 
SU(3) gauge theory is jllsu(2) /Msu(3) = 1-81(7). This result recalls analogous behavior seen in a 
different study of SU(2) and SU(3) vacuum in a constant external chromomagnetic background 
field [JTJ]. In Ref. [^|] numerical evidence that the deconfinement temperature for SU(2) and 
SU(3) gauge systems in a constant Abelian chromomagnetic field decreases when the strength of 



the applied field increases was given. Moreover, as discussed in Refs. [27, 41, 42], above a criti- 
cal strength y 'gH c of the chromomagnetic external background field the deconfined phase extends 



to very low temperatures. It was found [ ]41[ ] that the ratio between the critical field strengths for 
SU(2) and SU(3) gauge theories is \/ gH c I SU(2) /VgH c I su(3) = 2.03(17), in remarkable agreement 
with the ratio between the penetration lengths for SU(2) and SU(3). As stressed in the Conclu- 
sions of Ref. [r4l]|, the peculiar dependence of the deconfinement temperature on the strength of 
the Abelian chromomagnetic field gH could be naturally explained if the vacuum behaved as a 
disordered chromomagnetic condensate which confines color charges due both to the presence of a 
mass gap and the absence of color long range order, such as in the Feynman picture for Yang-Mills 



theory in (2+1) dimensions [43]. The circumstance that ratio between the SU(2) and SU(3) pene- 
tration lengths agrees within errors with the above discussed ratio of the critical chromomagnetic 
fields, suggests us that the Feynman picture of the Yang-Mills vacuum could be a useful guide to 
understand the dynamics of color confinement. 

References 

[1] M. Bander, Phys. Rept. 75, 205 (1981). 

[2] J. Greensite, Prog. Part. Nucl. Phys. 51, 1 (2003). 

[3] M. Fukugita and T. Niuya, Phys. Lett. B132, 374 (1983). 

[4] J.E. Kiskis and K. Sparks, Phys. Rev. D30, 1326 (1984). 

[5] J.W. Flower and S.W. Otto, Phys. Lett. B160, 128 (1985). 

[6] J. Wosiek and R.W. Haymaker, Phys. Rev. D36, 3297 (1987). 

[7] A. Di Giacomo et al, Phys. Lett. B236, 199 (1990). 

[8] A. Di Giacomo et al, Nucl. Phys. B347, 441 (1990). 



6 



Contribution title 



Alessandro Papa 



[9] V. Singh et al, Phys. Lett. B306, 1 15 (1993). 

[10] P. Cea and L. Cosmai, Nucl. Phys. Proc. Suppl. 30, 572 (1993). 

[11] Y. Matsubara et al, Nucl. Phys. Proc. Suppl. 34, 176 (1994). 

[12] P. Cea and L. Cosmai, Nuovo Cim. A107, 541 (1994). 

[13] P. Cea and L. Cosmai, Nucl. Phys. Proc. Suppl. 34, 219 (1994). 

[14] P. Cea and L. Cosmai, Phys. Lett. B349, 343 (1995). 

[15] P. Cea and L. Cosmai, Nucl. Phys. Proc. Suppl. 42, 225 (1995). 

[16] P. Cea and L. Cosmai, Phys. Rev. D52, 5152 (1995). 

[17] G.S. Bali etal, Phys. Rev. D51, 5165 (1995). 

[18] R.W. Haymaker and T. Matsuki, Phys. Rev. D75, 014501 (2007). 

[19] A. D'Alessandro et al, Nucl. Phys. B774, 168 (2007). 

[20] G. 't Hooft, in High Energy Physics, EPS International Conference, Palermo, 1975, edited by 
A. Zichichi (1975). 

[21] S. Mandelstam, Phys. Rept. 23, 245 (1976). 

[22] G. Ripka, hep-ph/03 10102. 

[23] A.A. Abrikosov, Soviet Physics JETP 5, 1 174 (1957). 

[24] H. Shiba and T. Suzuki, Phys. Lett. B351, 519 (1995). 

[25] N. Arasaki et al, Phys. Lett. B395, 275 (1997). 

[26] P. Cea and L. Cosmai, Phys. Rev. D62, 094510 (2000). 

[27] P. Cea and L. Cosmai, JHEP 11, 064 (2001). 

[28] A. Di Giacomo et al, Phys. Rev. D61, 034503 (2000). 

[29] A. Di Giacomo et al, Phys. Rev. D61, 034504 (2000). 

[30] J.M. Carmona et al, Phys. Rev. D64, 1 14507 (2001). 

[31] P. Cea et al, JHEP 02, 018 (2004). 

[32] A. D'Alessandro et al, Phys. Rev. D81, 094501 (2010). 

[33] G. 't Hooft, hep-th/0408183. 

[34] D.S. Kuzmenko and Y.A. Simonov, Phys. Lett. B494, 81 (2000). 

[35] A. Di Giacomo et al, Phys. Rept. 372, 319 (2002). 

[36] M.S. Cardaci et al, Phys. Rev. D83, 014502 (201 1). 

[37] N. Cabibbo and E. Marinari, Phys. Lett. B119, 387 (1982). 

[38] M. Campostrini et al, Phys. Lett. B225, 403 (1989). 

[39] R.G. Edwards et al, Nucl. Phys. B517, 377 (1998). 

[40] P. Bicudo et al, PoS LATTICE2010, 268 (2010). 

[41] P. Cea and L. Cosmai, JHEP 08, 079 (2005). 

[42] P. Cea et al, JHEP 12, 097 (2007). 

[43] R.P Feynman, Nucl. Phys. B188, 479 (1981). 



7 



